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DETERMINATION OF*THE STRESSES PRODU-OED- BY .THE LANDING 
IMPACT IN THE IKJLKHEADS OI' A SBAPLA-NB BOTTOM* 
By T. M. Darevsky , 

INTRODUCTION . 



The present report deals with the determination of 
the impact stresses in the "bulkhead floors of a seaplane 
"bottom. The dynamic problem is solved on the assumption 
of a certain elastic system, the floor heing assumed as 
a weightless elastic 1)63111 with c oncent rat ed • ma s s e s at the 
ends (due to the mass of the float) and with a spring 
which replaces the elastic action of the keel in the center, 
The distrihuted load on the floor is that due to- the hydro- 
dynamic force acting over a certain portion of the hottom. 
The pressure di str ihut i on over the width of the float is 
assumed to follow the Wagner law. The formulas given for 
the maximum "bending moment are derived on the assumption 
that the keel is relatively elastic, in which case it can 
he shown that at each instant of time the maximum hending 
moment i s • at the point of juncture of the floor with the 
keel. The "bending moment at this point is a function, of 
the half width of the wetted, surface c and reaches its 
maximum value when c is approximately equal to b/2 
where "b is the half width of the float. In general, 
however, for computing the tending moment the values of 
the "bending moment at the keel for "certain values of c 
are determined and a curve is draw.n. The illustrative 
sample computation gave for the stresses a result approx- 
imately equal to that o"btained hy the convent i onal .fact ory 
comput at i on . 



METHOD Of COMPUT AT I'ON 



Of the possible cases of the landing of a seapiane 
the most critical case will be that fo"r which the hydro- 
dynamic impact is most nearly .central (landing on the step). 

♦Report- -No. 449, of the- Cgntr'al -A-er o—Hydr odynami cal Insti- 
tute, Mosccsw, 1939". , •■ - 
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This case, therefore, will "be the Taasis of this discussion. 
In the cases of landing on the how or stern the length of 
the impact surfad'e'ls sufficient for-'the im;^act phenpmenon 
to he schematized as that of an infinitely long wedge on 
water. (Tlis prohlem'has heen solved hy Wagner ,( refer ence 
l)«) In the case of landing on the step, if the landing 
is effected on a calm surface, the length of the impact 
surface is small, and account must he taken of the fact 
that the float enters the water at a certain angle. In 
the most severe case, however - namely, landing on a 
large wave - the length of the impact surface exceeds the 
width hy at least 1.5 times. Hence, for this case the 
preceding schemati nation 1 s permi ssi-hio . - , 

In landing on the step when the ratio of the length 
of the impact' surface to the width is near 1.5, a correc- 
tion for finiten'ess of the surface is nevertheless desir- 
able. This correction may he made with the aid of the 
exper ime-ntal formulas obtained hy Pahst (reference 2) and 
Povitsky (reference 3). 

The number of main bulkheads ass-oclated vrith. the 
impact sur-face 'is not large. These bulkh.eads generally 
differ little from each other. For simplicity and with 
small error 'th6y will be assumed as' working under the 
same conditions. In accordance with the preceding assump- 
tion consider that the bottom of each of the bulkheads 
associated with the impact surface is loaded by a force 
P/n where P is the total impact force and n the number 
of bulkheads considered. In the case where certain of the 
bulkheads differ appreciably in stiffness it is sometimes 
possible to assume that the loading on the bulkhead bottoms 
is proportional to their stiffnesses. Since the bottom 
always has considerably less stiffness than the remaining 
structure (hull) of the float, consider the latter as abso- 
lutely rigid and only the bottom as elastic*. In this case 
in the landing on water the inertia forces will be trans- 
mitted mainly at the ends of the bulkhead floors. This 
pe rmits a considera t ion of the ma ss o f the entire float as 

*A mi sunder stant? ing may arise. On the one hand use is 
made of the Wagner solution which is applicable to a wedge 
of absolute stiffnes.s and,, on the other hand, the elasticity 
of the bottojn. The fact is that the hydr ody'namic pressure 
is determined by the. Wagner method. For a 7-bott.om float, 
however, the ehan-ge in-pressure for -s-mal 1- . change s in. the 
bottom 'shape due to elasticity is insignificant. .This 
change should be taken into account .only for practically 
flat bottoms. 



NAOA 'Technical Memorandum cITo . 1055 



3 



c one ent'ra"t'e'cL'"a't«"t-h'^e- ends of the floors*. If the mass of 
the entire float is denoted "by m, the mass associated 
vxith each end of the floor is m/2n. 

In the case of "eccentric impact, inst'ead of m, 

a = m i should Tie used where i is the radius of 

. .i.? t r^ ■ .. 

inertia of the^seaplane and- r "the.l^ever arm. relative to 
the center of gravity" of the seaplane. 

The centers of the floor 'beams will l)e acted oh lay 

the. elastic force of the "bottom with. t;he force transmitted 
mainly, through' the. keel. .This, force arises from the defor- 
mation of t.he TDott'op', on impact. The elastic force of the 
hot'tom is considered a.s the elastic force of a fixed keel . 
heam supported on intermediate elastic supports ( tha floors 
lying "beyond the impact surface). 

In order to take into account, to some extent, the 
action of the skin, the keel "beam and the bulkhead floors, 
with the a-ttached strips of skin to "be of a width of 40 
times the thickness, are considered together. 

Since the local action of the impact extends over a 
finite le-ngth of hottom, it will "be sufficient to consider 
the keel "beam to he fixed on three main "bulkheads to the 
left and the right of the impact surface.. Thus the part 
of th-e.'keel of interest may "be represented "by a "built-in 
"beam on four intermediate elastic supports (fig. 1). The 
rigidity of. these .supports k is determined "by the. formula 

n ! ^, whe"re / 2"b ' is the length of -the floar, - EI' t-he 

, -b3 - .-, -y . . , .$ 

sti.ff:nes-s (the- floors are considered as ^straigh^; "beams of 
constant stiffness). 

Th-e met/hod' of comput-ing "beams -on elastic suppoT-t s is 
desarijh'e-.d-.in de-.ta.il in the work of Umansky (refe.r,ence 4). 
By use ..of t-his n^e^thod-the maximum defleq-tion du-e ,to the 
conce;ntr,ated ;f orce at the -center of the earn P-.is-. 
determined' £ ■ 

^ : f- ' I 

*The "bull^he-ad fl oor s , assoclat ed with" the im.pact suVf ac.e . , 
only are 'considered. ' ' - " 
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K = K(BIg, K, AO . 

BIjj- stiffness of the keel 'beam 

K .stiffness of the supports (the ratio of the reaction 
of the supports to their deflection) 

Al spacing hetwee.n supports 

The elastic action exerted on the bulkhead floors of 
the v/etted surface hy the rest of the Ijottoin structure is 
thus replaced "by the action of an equivalent spring of 
stiffness K. Bach of the floors is acted on "by a spring 
of stiffness K/n. 

As a result the follouinri' scheme is s-rrived at. The 
elastic "bent heam (floor) v/ith the concentrated masses 
m/2a at the ends (due to the mass of the float) ha,s at the 
insto.nt t = 0 a vplocity Vq (initial velocity of immer- 
sion of float*). Thereafter the heaja is acted upon by the 
dis.tr ibuted load P/n (hj'-dr odynaaic pressure) and the con- 
centrated f or ce,s" . (elas t ic action of the bottom). At the 

center of the heam the fQrce p = where S is the 

n 

verticc.l displacement of the keel with respect to the chine 
and at the ends of the heara the forces are r/2 (fig. 2)** 

On account of the symmetry ono-half the heara 'liiay hd 
considered, the corresponding "boundary conditions being 
applied to the second half and the spring force taken to 
"be 5'/2 - that -is, the stiffness decre'ased "by one-half. 

The floors may he of the most varied shapes. Two 
•sharply differing shapes are shown in figures 3 and 4. 

• The elastic axis of -the floor is, ih general, curvedu- 
However, it ehall^he considered straight, since the' 
error made in thi s -^ s sumpt i on " shoiild be " insignificant ■ 
The elastic axis may be drawn as indicated in figures- 
5 and 4. Furthermore, the stiffness of the floor shall 

*!B'or landing on the step V^^. should be taken as the 
component of the landing velocity normal to the keel at 
the step. In the remaining cases msy be taken as 

the vertical eom'pbhent " of ' the landing velocity 
**The dotted line represents the stiff frame repTajsing ,the 
float hull. 



HACA Technical - Memorandum Ko. 1055 



"be considered as constant (mean, stiffness) althOMgh a 
variable stiffness offers no fundamental difficulties 
for the discussion. It aay "be noted that for modern 
floats there is a tendency to make the floor beam of 
almost constant height. In viev/ of the/ fact that the 
elastic axis is inclined at a certain angle to the 
surface of the water and the pressure of, the latter is 
not normal to the axis, the floor beam during the 
landing of the seaplane experiences both a transverse 
and a longitudinal impact. By decomposing the forces 
acting on the floor beam, including the inertia forces, 
into tuo components (along and perpendicular to the 
elastic axis) the transvers'e and longitudinal impacts 
may be separately ponsidered, the stress in the floor 
beam being determined as the sum of the. bending and 
compressive stresses. It 'is easy to show, hovrever , by 
siEiple computation that the compressive stresses will be 
negligibly snail by compar is on 'with the bending stresses 
and hence the longitudinal impact need not be considered. 
In exactly the same way, it is of no value to take into 
account the inclination ,of the elastic axis to the water 
surface, since this angle a will generally be 10 to 
20°, cos a = 0. 98430 — 0. 93869, and lience the corre- 
sponding components of the forces v;ill differ little from 
the values of^ the latter themselves, 

Shus the final scheme (figs, 5 and 6) will be some-» 
what simpler than that shown in' figure 2. J'lgur'e 5 shov;g 
the scheJie at the initial instant t = 0, figure 6 at the 
instant t v/hen the float is immersed someuhat in the 



wat er 
stiff 
S ince 
small 
r ight 



As in figure 2, the dotted line represents the 
frame replacing the uppev rigid part of the float* 
the mass of the beam (floor) shown m figure 5 is 

by comparison vith the mass concentrated at its 
end, the beam is considered as vreightless. 



The load distribution on the beam is determined vipb. 
the aid of the formula of Wagner (reference 5) for the 
hydrodynamic pressure on a submerged ii^edge: 



p(x,c) = 



Cl + V.)^ u 



1 + U - 



1 + - ^ 



r 



/I + u- 
1 + p. v • 



.3 



CU 



y 



du 
dc- 



1 
2 



u 



1 + u^ - 



-s-d + U2) 

_3 



- 1 



(1) 
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wher e 

p density of v;?ter 



G ha, If width of wetted surface of hottoni 
I longth of wetted surface 

and u IS deTierrained by the fornnila 
u = §f.an^. c / - iaiil] 



(2) 



P V angle 



Tl = T1 C ^ ) equation of oot^om profile descrilied rel-^tiv:- 
to the axes of coordinates shown m fiif^uro 7 

If the hottoji profile tiC^) nay he expres'sed m 
the form of the series* 

T1 (|) = Poe - Pll^ + ^zi'^ + . ■ . 

then u(c) is very simply determined as 

u(c) = ko^o + i^iP^c + ^'s^sc^ + 

wher e 

ir^ = C.S36, ki = 1.0, ks = 1.2V2, 

.13 - 1.5, ki = 1.6 96, Kg = i.c75 

and so forth** 

The linear loading on the oca-' tnercforo vill "b-^ 

o(x,c) = r-p(x,c) Al at [0,c] 
I C nt [c,o] 

where 

Al = 1 
n 

i length of impact surface 

h lengtii of h eoci ( ha If _v/ 1 d t h of fl ogt) 

*&ener?lly, for surfaces without discontinuities, tr.e 
profile IS sufficiently well represented hy an algehrpjc 
curve of the form 

^ = Poi - ^i-if 
**I?ho general formula for may he found m the A^rol. no 

Constructor's Reference Book, vol. II, 1958. 
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The equation of motion of the points of the elastic 
axis of the beam with linear loading a(x,c) io the 
following ^43, 

where a ie the density of the "beam. In this case the 
beam is weightless and therefore a = 0. On the other 

hand, the acceleration — ^ is a finite quantity owing 

dt 

to the finite mass at the end of the beam. Thus 

and the equation of motion assumes the form 

El^^,=q(.x,c) (c = c(t)). (3) 

with the boundary conditions 

1. y,(0.c) = 0, 2. Ery,"'(0,c)-k\y{b,c)-y(0,c)]=0. 
3. j;/(6.c) = 0, 4. E/y/'(b,c)-k[y(b,c)-y(0,c)]-^y,''{b.c)==0 , * 

where k = the stiffness of the spring 

Sn 

X. u« 1. .ad. Of t.. relation: ^ =g (| J+:|..^ , 

the fourth condition may be rewritten In the following 
form: m m 

A.£ly,'"(b,c)-k\y(b.c)-y{0. c)]-^^y"Ab, c) c'^- ^y'^b. c) c" =^0. 

The initial conditions are: 1. j> 0) = 0, 2. _y/ (^>, 0) c'j_5 = Vq. 

i>enote any function of x through p(x,c) by 

Pj(x,c), any function Pj^(x,c) by Pg(x,c), and so forth. 

From equation (3) there is obtained 

^(Pr (X, c) + C, (c) at [0. c], 
dx' \ PAc. c) + C,{c)at [c, b], 
il3L = if'? ^) + (c) X + (c) at [0, c], 

dx' I P, (c, c) + Ci (c) C+C, (c) + [Pi ic, c) + C, (c)\ (x - c) at [t,' *]. 
Pt(x,c)-[-^C,{c)x' + C,{c)x-{-C3 ic) at[0, c], 
P.ic.c)-}-^ Ci ic)c^ + C,ic) c + C, (c) + [P, (c, + Q (c) c+ 
+ Q (c)] (x_c)+[Pi ic c) + Ci at [0, b], 

*V(o account is taken of the weight of the concentrated mass 
since for a sufficient immersion of the float the mass Is 
compensated by the hydrostatic pressure and is moreover a 
small magnitude by comparison with the inertia forces. 



ay 

dx 
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Pi ix,c-{-^ C, (c) ;c» + ~ a {c) X' + C3 {c) x-^C, (c) at [0, c], 
P. (c, c)+-^C, (c) C +i- C, {c) c^- + C3 (0 c+Q IP3 (c, c) + 
+ i-C, (c) c' + C, {c) c + C, (c)] (X - c) + [P, (£:, c) + C, {c) c + 

With the aid of these relations and after the usual 
transformations the houndary conditions can he written as 
follows: 

1. Ct{c)=-P,{Q.c). 



C, {0 - ^ C, (0 = - EIP, (0. c) + Ub~ cf P, ic c) + 



+ -| (i, _ cy P, (c, c) - kbP, (0, c)-\-kib-c) P3 (c, c) - kPt (0, c)+ kP, {c. c). 

3. AC, (c) -f- C, (c) = -{b-c)P, (c, c) — P, (c, c). 

4. C/'(c)4-^ c"C/(c)-£/[Pac. ^) + Q(c)] + 
+ft[p«(^. c)-P4(0.c) + -^Q(c)c'+-^C,(c)c* + C3(0'^ + (^3(c, 0 + 



+ (c) 



j- P, (c, c) + C,(0 



(b — cf 



2n dc 



-2-^PAc.c) + PAc,c)-h^C," (c) c^-{-^C,"ic)c'+C,"{c)c + 



+ ^ Q' (c) + C ic) c 4- C (c) j + [ P, (c. c) - P, {c, c) + 
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From the first condition CsCe) is determined and 
from the second and third C^(e) and C,(c). Having 
determined Oi(c), C^Cc), and 03(0), the coefficient 
C4(c) is determined from the fourth "boundary condition 
and since this condition is a differential equation of the 
second order with respect to 04(0) the latter is obtained 
with two arbitrary constahts the ralue of which is found 
from two initial conditions. But 0^(c) need not be 
determined hecause the normal stress cr Is expressed 
through the bending moment in the expression for which 
04(0) does not enter. Thus by making use of the expree- 

sion for — ^ the following' formula is obtained for the 

bending moment M: 

f £/ [P, {X, c) + C, (c)x + C, (c)] at [0, c], 
\ EIlP,{c,c) + CAc)c + C,(c) + iPy{c,c)-\-C,ic))ix-c)]mt[c,b]. 
Trom the second ^nd third boundary conditions an 
expression is obtained for the values of Oj(c) and 

Ca(c) : 

C,(c) = ^- EIP, (0, c) + -g- kib-c) (c' -2bc - 2b') Pi {c, c) + 



-\-^kc(c- 2b) P, (c, c) - kbP, (0, c) + 

+ ft (6 - c) P, {c, c) - kP, (0, c) + ftp, {c, c)l ^ . 

J EI+^ 

C, (c) = ]^EIbP, (0, c)-{b- c) (eI- ^ + P, (c, c) - 



(5) 



^P,ic,c)-\-kb'P,iO,c)- 
—kb(p— c) P, (c, c) + kbP^ (0, c) - kbP^ (c. c) 1 



(6) 



ftA" 



Integrating successively with respect to x the 

function - p(x,c) results in the following expres- 
X I 

slons for P^Cx.c). Pa(x,c), P^Cxjc). and P4(x,c): 

'^-m% i[ ' [' +"'+-^ -Tine +"'=>]""M7fe?- 



I mill I ini^ll IIIIHIIII IHIHI^^HMHHIIIIBIIHI mil I 
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. . / , , , 3c* \ I* r2c'(i+tt»)+x« , /T~,~^ w, 
+0]/ i+^^--^j-r:fi^[ 3 V ^+«^--?-+ 

+ c(l+««)^arc sin-p^j— ^ |^_j,2_2<:«(l + «') + 
+ arc sln-^ + 3x.|/ - 

-*c-)=^ [Pf-^+^)x 

X \c {\-^u?)x(x^ + \-c'' (l + ««))arc sin -7=4== + 
L \. * / C|/ 1 -f- tt^ 



Substituting these expressions in formulas (6), (s), 
and (4) glTes the bending moment distribution along the 
beam. HoweTer, the value of the moment at the critical 
section is of interest. For a float bottom with weak 
keel (such as is usually the case in present day designs) 
the critical section will be at the point x » 0. 
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Thu« 

is found, where 

+ -g- (1 + J } arc Sin ^7^+ [ 1 + '^*+cu^c " 

-- 5r+]x^^+"^J 2(i+H^) 3-+— J- 

fc., . ftb' . kbc*-[ r94-10«» , 3 + 2tt» , dul 



+ - (2 + «*) J + c [ £/ — g- H g ^ 



■ «*)j+cj£/- 



+ |**cM4 + «')]ln(l-i^ 



(8) 



The value of u is, hoWeTer, sufficiently sBall 
(for a wedge with straight side's u > £ tan 3), so that 

TT 

in most of the terns of equations (?) and (8) u' and 

eu^ may be neglected, 
do 
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Finally, there is obtained 



o u(l + 

, 

720(1 + kt) 

24(1 + M,) L 



(9) 



A(c) = [ -[360TTEI + 10kc3(9Trb _ 8c) + 



ukc2(45TTh - 48o)] - ^ • xf61I[2"b(l - M-) + 5c(l + p.)] + 



khc[l0b^(l + M.) + 2lDc(l + 2|j,) + 0^(3 2ia.)]j + ^ (2EI + 

kbc^) In 2 + uc rgjji + kT3(T3^ + 

</ 1 + u - 1 24 



c 



)] In (l 1-^) 

\ 1 + u^/_ 



JJ (10) 



EIc + 

3 



The following should "be remarked. In the square 
■brackets of formula (lO), the second, third, and fourth 
terms, notwithstanding the presence of the factor u, 
may attain approximately a value of 20 percent of that 
of the first term. It may "be shown, however, that for 
the usual "bottom shapes and in the case of a weak keel, 
where the stiffness of the latter may he of the same 
order as the stiffness of the strongest floor beam, the 
algebraic sum of the preceding terms is small hy compari- 
son with the first term and constitutes not more than 

5 percent for c ^ — and not more than about 10 percent 

2 

for c > On this basis equation (lO) may be replaced 

by the equation 

A(c) = , — r- [SSOttSI + 10ko^(9TTb ~ 8c) + 

720(1 + M,) 

|j.kc®(45TTb - 48c)] ^ , (IC) 

EIc + 

o 

which will be sufficiently accurate for o <: i and will 

' 2 , 

give an error of the order of 10 percent for c > ^. 

2 
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By f-or-mulas 9.-) an.d (lO), -or (9.) .and (lO^) the stresses 
in the floor as a function of ,c aye ■ det ermined. To oh- 
tain the maximum stress,, a curve ^o^-^- drawn, a few 

values of the "bending moment "being -determined "by (9), and 
the maximum value of Mq(c) attained. This value should 
"be obtained, in general, for a cerfaln intermediate value 
of c. Thus, since the impact forcie increases with increase 
in c, the tending moment MqCc) will likewise at first 
increase. On the other hand, since with Increasing c 
the force of the hydrodynamic impact is concentrated 
usually at the edge of the wetted surface, it should he 
expected that after a certain instant of time the value 
of M^(c) should "begin to decrease, particularly since 

with increasing c the reaction of the keel increases 
(fig. 6). 

For V angles exceeding 15° the right side of equation 
( 1 ), and therefore also of equation (9), should "be multi- 
plied "by the correction factor u: 

u = 1 - JLi - 0.15 ^ - ^ In i 

where is the angle, in radians, of inclination of 

the bottom at the edge of the impact surface. 

Moreover equation (l) was o"btained for the conditions 
of the tw^~dimensi onal problem; that is, the wedge was 
assumed as infinitely long. Practically, equation (l) 
may "be used when the length of the wedge is Sufficiently 
large "by comparison with the width, at least one and one- 
half times as large. As has already "been remarked at the 
"beginning of this paper, for landings on the "bow and on 
the stern the length of the impact surface .is sufficiently 
great and may "be o"btained hy considering the landing on a 
calm surface. 

For landing on the step the most unfavorable condi- 
tion should he considered - namely, when the float lands 
on a wave. As experiments have shown, the length of the 
impact surface then exceeds the width "by 1.5 and more 
times. Since, for a ratio of length of wedge to width 
near 1.5 the error in computing the impact force on the 
assumption of two-dimensional flov/ is appreciable, a 
correction for finiteness is desirable. 



This correction may be made with the aid of the 
experimental formulas obtained by Pabst for the associated 
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mass 'o'f flat plates and' by Po'vi.tsky far 'the associated 

mass of wedges. For X >1.5, and usual "V angles, the 

'b 

formula 'of Pahst differs little from that of POvitsKy. 
B.oth 'of them will "be given. 

Let 'b 'be the half width of .the plate and I its 
length. The associated ma,ss M 'for the plate for the 
assirmption of two-dimensional flow is 



(11) 



the ' associated mass for the finite plate >l^is: 



M = ILp( Ic^ - c=^) 
2 

For an element of an infinitely long wedge the 
associated ma ssis 



(12) 



M = 



5— pic 

cos^3 



(13) 



where 

I length of element of wedi°:e 

c half width of wetted surface of wedge 

P V angle 



5 = i sin 2p 
3 



- 1 



(r denotes the Gam.iia function ) 

The value' of 6 also .may "be o'btained from the curve of 
figure 8. 



For the v/edge of finite length' l( 



/ -L ::^-l .425 , where 

Vsd" 

'd is the length of the wetted part of -the side of the 
wedge^ 

M = ^pic^ fi ^ ^\ ■ (14) 

cos^e V ly 
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The impact force for the wedge is exprfessed hy the 

formula - ---- - 




If the impact force, for the case of the two-dimensional 
prohlem is denoted by j^^^f true impact force for 

the finite wedge by P. , then 



dM^ 




= qp( c) 



dc 



Without taking into account the pressure distrihu- 
tion along the width of the wedge (the pressure distribu- 
tion is not identical for the infinite and finite wedges) 
a correction is made for finitenoss by multiplying the 
pressure determined by formula (l) by cp(c). In formula 
(9) it is then likewise necessary to multiply the right 
side by cp(c). By use of formulas (ll) and (l?.), there 
is obtained 



(15) 



1 - 



Sm + i. 
•npc^ c 



Prom (13) and (l4) is obtained 



1 - 



Sp c ^ c 



Thus, for comput ing. the bending moment M and the 
correction coefficient cp(c) for landing on the step, it 
remains merely to determine the length of the impact sur- 
face K It s'ho'uld be noted that, during the immersion 
of the float in the water, the length I is, in general, 
variable. It is nevertheless possible' to assume a certain 
constant value for I since,' in the most unfavorable case - 
namely, landing on a wave - the float almost instantaneously 
comes in contact with the wave over a considerable length 
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of surface which afterward increases only slightly. !For 
the same reason it is of no significance to make any 
correction for finiteness at the first instants of im- 
mersion. The correction should "be made when a sufficient 
width of the surface has already been immersed. 

The character of the sea surface permissihle for the 
landing of seaplanes varies with the. tonnage. The folloxir- 
ing mean values for the limiting height of the wave h 
in relation to aircraft of various tonnage ranges may "be 
given (reference S): 

First group G = 700 - 1500 kilograms, h = 0.15 - 0,3 motor 
Second group G = 3500 - 5000 kilograms, h = 0.8 - 1.0 meter 
Third group G = 8000 - 15000 kilograms,- h = 1.5 _ 1.8 meters 
Fourth group G = 20,000 and above, h = 2.0 meters 

With improvement in the seaway the preceding values 
may he raised somewhat. The length of the wave X is 
determined as a function of h by the curve of Zimmerman 
(reference 7). The xirave is sufficiently i^ell rep:^esented 
by a trochoid. The equation of the trochoid in parametric 
form is expressed by 

X = -A- CP + h. sin CO y = 2i — il cos CP 
2tt^ 2 22 

where X is the length of the v/ave and h its height. 
The trochoid also may be obtained by a simple geometric 
construction (reference S). Tho latter is indicated by 
the thin lines of figure 9. ?or the trim angle the mean 
angle of inclination of the tangent to tho wavo increased 
by 3° to 4° may be taken. For aircraft of not too great 
tonnage (3500 to 500C kg) the condition of the sea may be 
taken as a "calm surface" (the wave intensity is 2 points) 
which corresponds to an average force of the wing equal 
approximately to 3 meters per second (2.8 m/sec at a height 
of 2 m from ground level and 3.6 m/soc at a height of 12 m 
from ground level). Under these conditions the height of 
the wave is 1 meter and its length 11 meters. The angle 
of inclination of the tangent to such a wave varies from 
zero to 15 so that the trim angle may bo taken as 10°. 

To determine tho length of the wetted surface, draw 
the contour of the float with the chosen trim angle on 
the profile of the v/avo, immersing the step to the chine 
■(fig. 10). 'It is necessary, of course, to see that the 
tail of the float does not dig into the wave (this refers 
especially to floats with long tails). The maximum 
length of the wetted surfaoe is taken as the value of I. 
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As tosts have shown, in landing on the step the pres- 
sures at the loft end of the wetted surface are found to 
he nogllgihle. This may he explained hy the fact that 
hecause of the longitudinal curvature of the hottom the 
points sufficiently far removed from the step come in con- 
tact with the wfter at later instants of time. Moreover, 
the distance of these points from the center of gravity 
of the aircraft apparently has some effect, I'or this 
reason, for the purpose of greater accuracy, it is neces- 
sary to take for I not the entire length of the impact 
surface hut a smaller value - namely, the length of the 
part at which the hottom is almost cylindrical. 

If the width of the float is not too groat, it may 

ho assumed that the maximum impact force will occur at 
immersion of the chine. Since, however, the ratio of tho 
width of the wetted surfF.co to the length after the 
instant when the float is already sufficiently immersed 
in tho water varies little, the vr-luo of cp(c) may he 
assumed constant with tho width of the float taken as 
tho width of the impact surface. 

ILLUSTSATI^ra EXAMPLE 



The computation is conducted for a symmetrical 
landing on the step. The fundamental data reauired for 

the computation are tho following: — 

& weight in flight (2500 k^) 

landing velocity (30.6 m/scc) 

b half width of float at step (75 cm) 

A sketch of part of keel with bulkheads lying heyond 
the impact surface is shown in figure 11 (tho length of 
the wetted surface must ho initially determined). 

Ijj. moan moment of inertia of part of keel hoam under 
consideration (560 cm*) 

mean moment of inertia of floors lying directly 
beyond impact surface (90 cm"*) 

I mean moment of inertia of floors associated with 
impact surface (313 cm'^ ) 

W moan resistance moment of floors associated with 
impact surface (20.85 cn^ ) 

E elasticity modulus (7.2 X 10^ kg/cm^) 
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Computation 
Computation of K.- _ 48£/^ 48-7,2. 10» -90 

■ /» 140» — ^ 

Gomputatlon of. k.— For the deflection and angular 

displacement at section X of a "beam on elastic supports 
under the action of a force P applied at the center df 
heam, the' following formulas are obtained: (Reference 4 
makes use of the same notation.) 

f. = foA; - M,BJ - QoC; - y,DJ + N,, 

where , cp£> , Mq , and 0,^ are the deflection, 
angular dlsplaceiQeiit , hending moment, and shearing force 
at the origin of coordinates. 

The values A , B . C. D , N . A ' , B » , C • , D • , N • at the 
nth support (from the origin) are determined "by the recur- 
rence formulas: 

n-l ^ 
V 



1=1 

_A A 

— 2i*i^i^n-l> 

/=1 

A A 

1=1 

\ «~ 2EI ' «-'~ &E/ J 

n—i 
n-l 

1=1 

n-l 

/=1 
n-l 

^n' = 2M'Vl. 
1=1 



and 
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Place the origin of coordinates at the left built-in 
end. (See eketeh of part of keel under consideration.) 
Then yo='Po = 0 

and in the formulas for and <p^ there only eater 

the values C, D, H and B • , 0», N«. 

The values of and are obtained by equating 

to zero the expression for y and qp at the right end. 
Since y and qp for the right end are no other than 
^Bt Cp6» the magnitudes ------- 

»4. ^5. ^^1. ^a. 



^11 Cg. 



'3> 



4» Dj, J>s, D3, 



Da. -U.. H,. Ua. N,. N4, Ng are determined In suecesslon. 



Then 



40^ 



2.40,3-10' 



= 1,98-10- 



C, = ^,-xQ^,=^^^,--15,88.10'.l,98.10-^g^^=7.85.1^^ 



A 



A 



Cs = 5, - >cQ/3_, - xC,/3_, = 27^no^ - 1 5.88 . lOM ,98 . 1 0-« . q,^^\q-, ■ 



15,88- 10'. 7,85- 10-" 



270' 



6 -40,3 -10' 



= 12,05-10 



A A A A 
C4 = xCj/4_i *CJi_^ xCg / 4-3 = 



390« 



2-40,3.10' 



"VW 310' 
-15,88-10'.l,98.10-«-^^ggi^-15,88.10'-7,85.10-". e.^0^^ 



— 15,88- 10». 12,05. 10-" 



40' 



6 -40,3 -10' 



= — 20,093-10-" 



A A A A A 



430« 



2-40,3-10' 



— 15,88-10'.l,98-10-« 



390" 



6.40,3-10' 



Vint ons 

15.88- 10'-7,85. 10- - 15,88. 10»- 12,05- 10- ^^^g^-f- 

-f 15,88 • 10* - 20,093 - 10-' = - 73,3 - lO-. 

In the same manner, by the preceding formulas the 
remaining magnitudes are computed. The results are given 
In the following table: 





1 


2 


3 


4 


5 


c 




1.98- 10-* 


7,85-10-® 


12,05-10-® 


— 20.093 -I0-* 


—73,3. 10-® 


D 




26,45- 10-® 


210,8-10-® 


14510-10-® 


19570-10-® 


20350.10-* 


N 




0 


0 


1016-10-®? 


2168-10-®? 


3682-10-®.P 


B' 












-1,504-10-® 


C 












—73,6.10-® 


N' 




- 








37,7.10-'-P 
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Now Mq, are determined from the equations: 

y-, 73,3 . lO-'M, — 20350 • 1 Q-'^Q,, -f 3682 • IQ-" P, 
9,, = 1 ,504 • lO-^Afo + 73,6 • 10~^Q„ + 37,7 ■ lO"' P. 



There Is obtained 



Qo = 



3682 —20350 
37,7 73,6 



73,3 
1,504 

73,3 
1.504 



20350 
73,6 



73,3 
1,504 



3682 




37,7 




20350 




73,6 





- 29,2 P, 



A> = 0,0769 P. 



Then the maximum deflection is found. The coordinate 
of the center of the beam is 



A:p = 215 cm 



Since 



A A 



= 2-W?W-^^'««-10^-^-9«-^«"''6TiflO^- 
- 15.88. 10-7,85. 10- = 37,7. 10 -6. 



6. 40,3. 10'' 



A A 

2158 



1758 



g:^0;3-jQ, -15,88.10=. 26,45.10-^:^,. 

135^ 



15,88. 10*. 210,8- lO-**. 



6- 40,3- 10 

:0 



. = 3682.10 ^ 



and therefore = 29,2 -P. 37,7. 10 — 0,0769. P- 3682. 10-'' = 845. 10 « P. 

K is determined A'= -i-. 10" = 11,83.10" . 

845 cm ' cm 

Since five bulkheads are associated with the wetted 

surface a' i 

k = -^ = 1,183-10- ifi . 
10 ' em 



Determination of the function u(c). - To determine 
the function u(c), it is necessary first to write down an 
analytical expression for the bottom profile. In this 
case the bottom profile consists of a section of a straight 
line and part of a circle, but it is more convenient to 
express it by a curve of the form Ti = - 3n_il'^- 

For n = 4, this curve very accurately represents the 
profile. The coefficient Pq and 3^ are determined 

by two points of the profile. 
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Denoting by 1^ and i^, respectively, the coordinates 

of the end of the straight- line and the .end of the curvi- 
linear parts of the profile gives from the sketch 

£j = 41 cm,,., -/ji = 14,4 cm , 
Sj = 75cm, ifi, =:20 cn. 

whence ^^^-^^ ^ 14.4.75«-20.4H ^ 

41. 75*- 75. 41* 

p. T.,£,-r, ^, 20.41-14.4-75 p 03.5 ip-n [JLl 

M7^«2V ~ 41.75*-75.41« -"•^'^^^ |„3| 

and since A, — 0.636, fci= 1,500, 

therefore' 

a(c)=-A!jp„ — /fe,,3sc' =0,636. 0.37 --l,5-0,2375.10-«<;»= * 
" . • •. =0,235— ;0,356. 10-" r'. 

Computation of I, Al , M-Cc), ^o-- To determine I, 

the wave is- -drawn to ^ certain scale in the form of a 
trochoid of height of 1 meter and length ll meters. To 
the same scale the contour of the float is drawn and 
superposed oh the wave with a trim angle equal to 10® 
(fig. 10). Then hy the same procedure as indicated 
previously, there is ohtalned from the sketch 

' ■ / = 200 cm 

*/ I 200 ^_ 
A/ = — = —-=40 CO 
n 5 

' :. , " ' ^ , , 0,001028 kg 

Assume for p the value p = T» 

^ g cm** 

There is obtained tx = p^ c' = 0,001028 . M|^c«= 1,47- 10-«.c« 

V^o = »^„ocSin 10° = 30,6.0.17365.jfj^5,3^. 
/' 2200 , 

where _ m, = m ^ _^ ^, s= ~j- Kg 

- Collection of data required for the computation. - 

f» = 75 cm. A/ =40 cm~ 

0,001028 kg 1 J 

p = — i- a- » ti = 1,47. lO^.c* • • 



u= 0,235 = 0,356- lO'^c', £"/=:22,57- lO"-^-, 
A =1.183.10'^. V/„ = 5,3-rfT. 



cm* 



Determination of A(c) and H(c). - By formulas (9) 
and (-10) are computed the values of M^^^^ ^ number 

of values of c. A curve is then drawn from-which the 
maximum value Mq(c) Is determined. 
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Por c = ^ = 37,5 centimeters, there is obtained 

IJ. = 1,47x10"^ X_37.5^ = 0.2065 

u = 0,2-35 - 0.356 X 10~^ X 37.5^ = 0.2162 

A (37.5) = -2.56 

Hq(37.5) = 8.2 X lO'* kilogram-centimeter 

By computing also IIq(c) for the points c = 20, 30, 

50, 75 centimeters, a curve ii = 1Iq(c) is ohtained (fig, 12), 

Jrom the figure it is seen that Mq(c) reaches a maximum 

approximately for c = — and Mq(c) = 8.2 X 10** kilogram— 

cent imet er . 



1 




/ i 


;:::::=Te-| 


2 X 2200 


+ 200 1 


3.14 X 1.028 X 10"^ X 


75=^ 75 1 


v 





ne.±.£l.!Biaaiiaa_fif._t.il£L_s.ir as.aaa a.— Since 9 X maximum 

Mq(c) = 6.28 X 10* kilogram-centimeter and the resistance 
moment of the floors V = 20,85 cuoic centimeters, there 
is ootained 

cp X max Ho(c) ^. . kg 
0 = = oU , 1 — 

¥ mm^ 

It may he noted that the present— day factory compu- 
tation (the so— called ''static computation") of the floor 
as an element of the transverse system gives tJ = 33 
kilograms per square millimeter. 



Translation by S. Heiss, 
National Advisory Committee 
f or A er onaut i c s . 
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Jigs. "1,2,3,4,5 



Impact surface 



Floor of ■bulkhead 



JPigure 1. 




Figure 2. 




Figure 5. 
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Pigs. 6,7,8,9 




Jigure 9. 
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figs. 10,11,12 
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Figure 10, 
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§.1. §-4 Wetted s-urface L B-4- 

Pigure 11, 




Figure 12. 



